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ABSTRACT 

Planets can emit polarized thermal radiation, just like brown dwarfs. We present calculated 
thermal polarization signals from hot cxoplanets, using an advanced radiative transfer code that 
fully includes all orders of scattering by gaseous molecules and cloud particles. The code spatially 
resolves the disk of the planet, allowing simulations for horizontally inhomogeneous planets. Our 
results show that the degree of linear polarization, P, of an exoplanet's thermal radiation is 
expected to be highest near the planet's limb and that this P depends on the temperature and 
its gradient, the scattering properties and the distribution of the cloud particles. Integrated over 
the disk of a spherically symmetric planet, P of the thermal radiation equals zero. However, for 
planets that appear spherically asymmetric, e.g. due to flattening, cloud bands or spots in their 
atmosphere, differences in their day and night sides, and/or obscuring rings, P is often larger than 
0.1%, in favorable cases even reaching several percent at near-infrared wavelengths. Detection of 
thermal polarization signals can give access to planetary parameters that are otherwise hard to 
obtain: it immediately confirms the presence of clouds, and P can then constrain atmospheric 
inhomogeneities and the flattening due to the planet's rotation rate. For zonally symmetric 
planets, the angle of polarization will yield the components of the planet's spin axis normal to 
the line-of-sight. Finally, our simulations show that P is generally more sensitive to variability 
in a cloudy planet's atmosphere than the thermal flux is, and could hence better reveal certain 
dynamical processes. 
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Studying the thermal emission of exoplanets 
has only recently become possible with the di- 
rect detections of young gia nt planets in wide 
orbits around their st ar (e.g. iMarois et al. 20081: 
Lafreniere et al. 20081) and secondary eclipse de- 
tecti ons of transiting exopl anets in tight orbits 
(see Deming fc Seager 2009t ). These hot planets 
emit most of their radiation at near-infrared wave- 
lengths and flux measurements at different wave- 
lengths are used to constrain properties of their 
atmospheres. However, little attention has been 
given to the information contained within the po- 
larization signals of this planetary thermal radia- 



tion. 

Incident starlight will usually get polarized 
when it is scattered by gases and particles in 
a planetary atmosphere. When integrated over 
the planetary disk, the reflected starlight can 
yield a significant net degree of pola rization (e.g. 
Seager et all 120001 IStam et al.ll2004l ). Also polar- 
ization in stellar atmospheres is well-known (e.g. 
Harrington 1970) and a dozen po larized brown 
dwar fs have been identified (e.g. Menard et al.l 



2002). Thermal radiation that is emitted by a 
body can get polarized upon scattering. For a net 
polarized thermal signal of a spatially unresolved 
body such as a brown dwarf or an exoplanet to 
be observable, not only scattering particles are 
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required, but also an asymmetry in the body's 
disk. Otherwise, the polarized thermal signals 
from different parts of the disk will cancel each 
other completely. With an asymmetric disk, the 
cancellation will be incomplete and a net polar- 
ization signal remains. 

For the polarized brown dwarfs, a plausi- 
ble source of asymmetry is the flattening of 
the body due to it s rota tion, as advocated by 
Sengupta fc Marlev ( 201dt ). Very recently, these 
authors have extended their work on polar- 
i zation of flattened objec ts to planets as well 
(|Marlev fc SenguptalEoil . In a planetary atmo- 
sphere, sources of asymmetry could arise from hor- 
izontal inhomogeneities in temperature or cloud 
thickness. For instance, strong zonal winds can 
cause banded structures, as found on all giant 
planets in our solar system. Furthermore, rings 
that obscure or shadow part of the planetary disk 
would cause asymmetries. 

Here, we will use simulated signals of horizon- 
tally inhomogeneous planets to present processes 
that polarize thermal planetary radiation, to ex- 
plore parameters that determine the strength of 
these thermal polarization signals, and to discuss 
the value of infrared polarimetry for the charac- 
terization of exoplanet atmospheres. 

2. Our numerical model 

Light is fully described by a flux vector F = 
(F,Q,U,V), with F the total flux, Q and U the 
linearly pola rized fluxes, and V the circularly po- 
larized flux (Han sen fc Travis! Il974ft . Parameters 
Q and U are defined with respect to a given refer- 
ence plane. For atmospheric calculations we have 
a local reference plane with axes parallel and nor- 
mal to the local planetary horizon. In the follow- 
ing, we will neglect V, and express the degree of 
polarization, P, using 



P 



VQ 2 + u 2 



(1) 



The angle of polarization, %, wi th respect to the 
refere nce plane is defined as (see lHansen fc Travis 
19741) 

tan2x = L7/Q. (2) 

We calculate the radiative transfer in a locally 
plane-parallel, vertically inhomogeneous planetary 



atmosphere for a range of emission angles usin g 
a doubling-adding method (jWauben et al. l ll994lV 
which fully includes all orders of scattering and 
polarization. Each atmosphere consists of 40 lay- 
ers, equally spaced in log pressure, with pres- 
sures p between 10 -6 bar (top) and 5 bar (bot- 
tom). In order to make comparisons between dif- 
ferent model atmospheres we use ad hoc temper- 
ature profiles, assuming hydrostatic equilibrium, 
which are uniform in temperature T for p < 10~ 3 
bar and p > 1 bar. In between these pressures 
the near-infrared emission originates and we as- 
sume dT/dhip is constant there, pivoting around 
T = 1500 K at p — 33 mbar for different tem- 
perature profiles. We perform our calculations at 
Ai = 1.05 /im, which is in the continuum (Y- 
band), and at A2 = 1.11 fini, in a water vapor 
absorption band. The water absorption optical 
depths ar e calculated using th e HITEMP 2010 
database ( Rothman et al. 2010l) . assuming Voigt 
line shapes and a volume mixing ratio of 5-10~ 4 . 

At these infrared wavelengths, the scattering 
optical thickness of the gas molecules is negli- 
gible, but radiation can be scattered by larger 
cloud/dust particles. We calculate the scatter- 
ing matrices of the particles in our model at- 
mospheres using Mie-sca t tering as described by 
de Rooii fc van der Stapl ( 19841 ) (thus assuming 
spherical particles). The particles that we use are 
small compared to the wavelength and have unity 
single scattering albedo. The scattering of these 
particles can be described as Rayleigh-scattering, 
with a nearly isotropic flux phase function and a 
bell-shaped polarization phase function reaching 
almost 100% polarization at a single scattering 
angle of 90°. The particles are assumed to be in 
one of the atmospheric layers and we assume iden- 
tical cloud optical thicknesses at Ai and A2, which 
only have a difference in wavelength of 0.05 /xm. 
We also performed calculations using larger, non- 
Rayleigh-like, particles, which will be discussed in 
Section 3.4. 

3. Spatially resolved polarimetry 

To understand disk-integrated polarization sig- 
nals of exoplanets better, we first discuss polariza- 
tion signals that are spatially resolved across the 
disk. 

Generally, singly scattered radiation is unpo- 
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larized for scattering angles equal or close to 
0° (assuming unpolar ized incoming radiation, 
Hansen fc Travis! 1974 ). Hence, thermal radia- 
tion emitted by low atmospheric layers that is 
scattered upward by particles in upper layers will 
have a very low degree of polarization. Radiation 
that is emitted in a high altitude layer that it- 
self contains scattering particles can also be scat- 
tered in the upward direction over a 90° angle 
and be strongly polarized, in particular when the 
scattering is Rayleigh-like. However, because the 
scattering particles in the upper layer will receive 
emitted radiation from all directions equally, the 
net polarization of the scattered radiation that 
emerges from the layer will equal zero. Hence, at 
the center of a horizontally homogeneous plane- 
tary disk, P will usually be low, regardless of the 
temperature profile. 

At the limb of a planetary disk, thermal radi- 
ation that is emitted by low atmospheric layers 
will be scattered by particles in upper layers to- 
wards the observer over angles near 90°, yielding 
high values of P. Radiation that is emitted in 
high altitude layers that contain scattering parti- 
cles themselves will yield relatively low values of 
P, because of the large range of scattering angles. 
P will not be zero mainly because of contributions 
of radiation that has been scattered twice or more. 

The degree of polarization at the limb thus de- 
pends strongly on the temperature structure of the 
atmosphere. This is related to the limb darkening 
effect, which is illustrated in Figure 1. For a given 
optical path length, radiation arriving at an an- 
gle towards the scatterer will come from a higher 
layer in the atmosphere than the light coming from 
directly below. If there is a temperature gradi- 
ent, that difference in altitude will correspond to 
a difference in temperature and hence a differ- 
ence in thermal flux. If the temperature decreases 
with decreasing pressure, most thermal radiation 
emerging from the limb of the planet will have 
been emitted in the lower atmospheric layers, and 
will have been scattered at high altitudes towards 
the observer at angles close to 90°, where P of 
Rayleigh scattered radiation is highest. If the tem- 
perature increases with decreasing pressure, most 
of the thermal radiation that is observed at the 
planet's limb will have been emitted and scattered 
in the upper layers. In this case, the large range 
of scattering angles will yield relatively low values 



of P (see Figure 1). If the temperature gradient 
equals zero, the radiation field around the scatter- 
ers is more symmetric, even at the limb, and the 
net P of radiation emitted towards the observer 
will be close to zero. 

3.1. Dependence on temperature gradient 

Figure 2a shows the normalized polarized flux 
—Q/F for model atmospheres with different tem- 
perature profiles and a high altitude cloud layer 
with an optical thickness of 0.1 at Ai and A2 for 
small and large emission angles. In these one- 
dimensional model calculations, U — 0, and hence 
I Q/F I = P. A positive (negative) —Q/F indicates 
that the radiation is polarized parallel (perpen- 
dicularly) with respect to the local horizon. The 
curves in Figure 2a show the strong influence of 
the temperature gradient on P and that a positive 
temperature gradient changes the polarization an- 
gle X by 90°. 

It is also clear that for positive temperature gra- 
dients, P is highest for A2, the wavelength in the 
water vapor absorption band. This can again be 
understood in terms of the effect shown in Fig- 
ure 1. In the continuum, the radiation that arrives 
at the scattering particle originates deeper in the 
atmosphere than in the water band. Because of 
our constant temperature gradient, the tempera- 
ture difference between radiation coming towards 
the particle from directly below and that coming 
at an angle will be roughly similar for both wave- 
lengths. However, such a temperature difference 
corresponds to a relatively larger difference in flux 
for lower temperatures, meaning that radiation ar- 
riving from the colder atmospheric layers probed 
in the absorption band will have relatively more 
contributions from directly below and hence will 
give a larger P at the limb. For negative temper- 
ature gradients, —Q/F is negative, which shows 
that here the emerging signal is dominated by ra- 
diation arriving at an angle (see Figure lb) that is 
scattered multiple times. In the continuum, there 
is now relatively less radiation arriving from di- 
rectly below, giving less positive contributions to 
the net —Q/F value. Hence, P is higher in the 
continuum than in the absorption band for atmo- 
spheres with a negative temparature gradient and 
high clouds. 
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3.2. Dependence on cloud optical thick- 
ness 

An optically very thin cloud won't affect the 
outgoing radiation much as the cloud is al- 
most transparent to the outward-going radiation. 
Hence, not much radiation is scattered and P will 
be low. With increasing cloud optical thickness, 
P will increase, because more radiation is scat- 
tered as the transmission of the clouds decreases. 
On the other hand, the contribution of multiple 
scattered radiation will also increase with increas- 
ing cloud (scattering) optical thickness. This is 
because the free path length of the radiation can 
become less than the distance between scattering 
particles. Multiple scattered radiation arriving at 
the observer will have been scattered at a range of 
angles. At the limb, this means that the observed 
radiation has not only been scattered over angles 
that produce a high P, but it will also have been 
scattered over angles that produce low P. For 
small emission angles, such as observed near the 
center of the planetary disk, multiple scattering 
can actually increase P, because it adds radia- 
tion that has been scattered at angles where P is 
large. For very large optical thicknesses, virtually 
all the thermal radiation will originate from the 
cloud and be multiple scattered within the cloud 
itself and P will be relatively low. These various 
effects of the cloud optical thickness can be seen 
in Figure 2b: first P increases with cloud opti- 
cal thickness as more radiation is scattered, then 
P decreases as multiple scattering becomes im- 
portant, and finally P stays roughly constant at 
large cloud optical thicknesses as virtually all ra- 
diation originates from multiple scattering within 
the cloud itself. 

3.3. Dependence on cloud top height 

In Figure 2, the cloud is located high in the 
atmosphere. Lowering the cloud does not change 
P much until the gas absorption optical depth is 
comparable to the cloud scattering optical thick- 
ness. For optically thin clouds, P gradually van- 
ishes as the cloud descents and gaseous absorp- 
tion starts dominating over scattering. For thicker 
clouds, P at the limb will first rise slightly be- 
fore vanishing as multiple scattering effects are re- 
duced. Because gas optical thicknesses depend on 
wavelength, cloud top heights can have large ef- 



fects on P's spectral behavior. For instance, in 
gaseous absorption bands, P can be zero when 
clouds are below the altitudes where gaseous ab- 
sorption is high, whereas it can be non-zero at 
continuum wavelengths. Hence, a case with low 
clouds can leave peaks in the polarization spec- 
trum at certain wavelengths, whereas in a case 
with high clouds these wavelengths represent a lo- 
cal minimum in the polarization spectrum. 

3.4. Dependence on scattering properties 

In the previous calculations, we assumed par- 
ticles that were small relative to the wavelength 
(i.e. they are Rayleigh scattering). The polar- 
ization signal of a planet will, however, depend 
on the particles' scattering properties. With in- 
creasing particle size, the maximum value of the 
polarization phase functions usually decreases 
and polarization direction changes occur (e.g. 



Stam et al.l 120041 ) . Calculations with Saturn-like 
pa rticles, whose scatt e ring p roperties are derived 
bv iTomasko fc Doosel ( 1984h . and 1-^rni spherical 
Mg2SiC>4 particles show that P at the limb will 
decrease to less than 10% and that the sign of Q 
will be opposite to that shown in Figure 2. In 
both of these cases the single scattering albedo of 
the particles is close to unity. 

3.5. Three-dimensional effects 

The above calculations were performed using 
one-dimensional models, which are applicable to 
atmospheres that are locally horizontally homoge- 
neous. However, large local variations in e.g. tem- 
perature or composition can give rise to adjacency 
effects that affect a planet's polarization signal. 
For instance, a hot spot in an atmosphere could 
be surrounded by a ring of high polarization re- 
sulting from radiation that is emitted by the hot 
spot, which is subsequently scattered by the colder 
atmosphere around it. Such three-dimensional ef- 
fects can be the subject for later study. 

4. Disk-integrated polarimetry 

To simulate exoplanet signals, we integrate the 
spatially resolved fluxes across the planetary disk, 
using a grid with a 2° resolution in latitude and 
longitude. The reference plane for the polarized 
fluxes Q and U depends on the location on the 
planet (jWauben et al.lll994l) . Before integrating 
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these fluxes, we thus have to rotate the local 
flux vectors to the planet's global reference plan e 
(|Hovenier fc van der Melll983l : IStam et al.ll2006h . 
which we align with the planet's spin axis. For 
spherically symmetric planets, the rotation and 
subsequent integration yields a complete cancel- 
lation of both Q and U of the planetary thermal 
radiation (see Figure 3). 

We have performed calculations for five differ- 
ent asymmetric planets, all illustrated in Figure 3: 
spherical planets with a band, a spot, obscuring 
rings (without ring shadows), a day- night differ- 
ence, and a flattened (ellipsoidal) planet with 
a horizontally homogeneous atmosphere. Obvi- 
ously, the parameter space for these calculations 
is enormous. Here, we present a limited num- 
ber of cases to provide an indication of the range 
of polarization signals of inhomogeneous planets 
and to draw some qualitative conclusions. We 
model our horizontally inhomogeneous planets us- 
ing two different atmospheric profiles, and hence 
two different one-dimensional radiative transfer 
calculations (at several emission angles). Depend- 
ing on the case (band, spot, etc.) and viewing 
geometry, we assign one of the two profiles to 
the grid points on the planet. We use a tem- 
perature gradient of 300 K/lnp everywhere on 
the planet, but in the clear parts of the planet, 
temperatures are 250 K lower than in the cloudy 
parts, which i s a modest temperatur e contrast for 



hot Jupiters ([Showman et al. 2009). In best-fit 



model calculations of the directly i maged planets 



around HR 8799 and 2MI207 by Barman et al 



(|2011bl lah ;he temperature gradient at pressure 
higher than 0.1 bar is on average steeper than the 
300 K/lnp used here, whereas below this pressure 
the temperature gradient is more shallow. The 
tempe rature profile shown in iMarlev fe Sengupta 
(|201lh shows very similar behaviour for T e ff=1000 
K, 3=30 m s -2 and / se d=2, although temperature 
gradients are slightly less steep around 0.1 bar. 
Also best-fit atmospheres of transiting exoplanets 
have temperatures gradi ents around 300 K/lnjo 
in th e lower atmosphere (Madhusudhan fe Seager 
20091 iMadhusudhan et al.l l2011bh . In the near- 
infrared, most of the thermal emission originates 
from the region in the atmosphere between 0.1- 
1 bar for positive temperature gradients (e.g . 



not seem unreasonably steep. For the cloudy part 
of the atmosphere, a Rayleigh scattering cloud 
layer with unity scattering optical thickness is 
placed at the top of the atmosphere. 

Our calculations o f ellipsoidal planets can b e 
compared to those bv lSengupta fe Marlev ( 2010l ). 
For a homogeneous cloudy atmosphere and an in- 
clination angle of 90°, our P at Ai as a function 
of oblateness compa r es we ll with the /-band P of 
Sengupta fc Marlevl (|2010t ) for e.g. T cS = 1800 K 
with log(g) = 4.5 and / Sfi ri=2 (see t heir Fig- 
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Showman etalll2009l: IMadhusudhan et alj2011bh . 
and hence our assumed temperature gradient docs 



ure 1). Like ISengupta fc Marlevl (|2010l ). we can 
reach values of P of several percent for extreme 
oblatenesses at Ai, and roughly twice that value 
at A2. Longer wavelengths will give only slightly 
lower Ps, given identical temperatures and opti- 
cal depths, as is shown in Fig. |4] (the cloud optical 
thickness is kept constant with wavelength). The 
figure also shows the effect of low clouds, which 
results in inverted polarization spectra. Also note 
the large difference in polarization spectra be- 
tween the two cases, whereas the flux spectra are 
very similar. Calculations with l-/im Mg2Si04 
particles, with a cloud optical thickness of unity 
at 1 fim and wavelength-variations over the plot- 
ted wavelength range determined by Mie theory, 
result in very similar polarization spectra, scaled 
down by a factor of ^5. 

For the planet with the equatorial cloud band, 
we varied the band's latitudinal width and the 
planet's inclination angle. A clear planet with a 
dusty band shows a maximum P (0.5% at Ai and 
2% at A 2 ) for a 40°-60° wide band. For a given 
band width, P decreases with decreasing inclina- 
tion angle, to increase again slightly when incli- 
nation angles get larger and the band reaches the 
limb. When the planet is seen exactly pole-on, 
P = again, as expected. For a dusty planet with 
a clear band, P can reach 0.5% at Ai and 4% at 
A2, depending on the inclination angle. 

Except for very special geometries, rings will 
usually obscure part of a planet either directly or 
through shadowing. We simulated the presence of 
cold, optically thick rings for a range of inclination 
angles by modeling the obscured part of the disk 
with a 200 K blackbody (see Figure 3). For a ring- 
planet radius ratio similar to that of Saturn, P is 
at most 0.3% at Ai and 0.8% at A2 at a ring incli- 
nation angle of ^45°. If the rings extend further 
outwards, and hence cover a larger fraction of the 
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planet's disk, the maximum value of P increases 
by a few tenths of a percent. Note that we only 
modelled the obscuration here. The rings them- 
selves can scatter the thermal radiation from the 
planet and might also emit significant polarized 
radiation themselves if they are hot and optically 
thick. A three-dimensional model of the planet 
and rings would be needed to model the interplay 
between planet and rings. 

All three cases above have the planet's spin axis 
as the axis of symmetry. As a result, the disk- 
integrated polarized flux U equals zero and the 
angle of polarization \ is either parallel or perpen- 
dicular to the spin axis, depending on the parti- 
cles' microphysical properties and the temperature 
profile. We now briefly discuss two cases without 
this type of symmetry. 

We simulated an equatorial cloudy hot spot by 
a dusty square of 20° x 20° (latitude x longitude) 
(see Figure 3). As the planet rotates, F varies 
mildly, while P and x vary significantly (see Fig- 
ure 4). Here, P reaches 0.1% at Ai, and 0.6% at 
A 2 . 

The planet with the day-night difference has 
a dusty and a clear hemisphere (split along lon- 
gitude lines, parallel to the terminator). As the 
planet rotates, P depends on which parts of the 
hemispheres are in view. The maximum P equals 
0.6% at Ai and 4% at A2, and P — when only 
one hemisphere is in view, x varies similarly to 
the spot case (Figure 4). Again, flux F varies only 
mildy compared to P. 

5. Conclusions and opportunities for exo- 
planet characterization 

We have calculated polarization signals of ther- 
mal radiation that is scattered by particles in ex- 
oplanet atmospheres. Important parameters that 
determine the degree of linear polarization P are 
the particles' polarization phase function, the op- 
tical thickness and the altitude of the particles, 
and the temperature and its gradient. 

For spatially resolved planets, P is usually high- 
est near the planet's limb. For spatially unre- 
solved planets, inhomogeneities on a planet's disk 
can cause a net polarization signal. In the cases 
we have considered (i.e. a band, a spot, obscuring 
rings, a day-night difference and a horizontally ho- 
mogeneous flattened disk), P was typically above 



0.1%, and values of several percent were reached in 
favorable cases. Combining different effects, e.g. a 
band on a flattened planet, could increase P even 
further. With a high cloud and a positive temper- 
ature gradient, P is higher in gaseous absorption 
bands than in the surrounding continuum. Unfor- 
tunately, little flux is emitted in those bands in 
this case. Together with telluric absorption, this 
will make the water absorption bands of planets 
with a positive temperature gradient less suitable 
for detection of infrared polarization in exoplan- 
ets using ground-based telescopes. Fortunately, P 
is at most a factor of a few lower in the atmo- 
spheric windows, whereas the flux can be several 
magnitudes higher, giving more opportunities for 
detection there. 

Currently, exoplanetary thermal radiation is 
detected either during secondary eclipses of tran- 
siting planets or through direct imaging of planets 
at large orbital distances. In the former case the 
combined light of the star and planet is measured, 
and although these usually tidally locked planets 
are expected to display large temperature inho- 
mogoneities, measuring a 1% polarization signal of 
the planet will be very challenging, as sensitivities 
of ~10 -6 need to be reached. Indeed, direct imag- 
ing, where the planet is spatially resolved from its 
star, promises to be a more suitable method to 
detect polarized exoplanets, especially with new 
instruments like Gemini/GPI and VLT/SPHERE. 

There are several reasons why planets like those 
in the HR 8799 system might be more prone to 
producing polarized signals than brown dwarfs, 
which can be polarized by a few percent. Firstly, 
the surface gravity of planets is lower than that of 
brown dwarfs, making them more flattened for a 
given rotation rate (see also iMarlev fc Sengupta 
201ll ). Secondly, broadband flux measurements 



suggest the HR 8799 planets to be more cloudy 
than most brown dwarfs, wh ich might be a com- 
mon fe ature of young planets ( Madhusudhan et al.l 
2011a[ ). Thirdly, the effective temperatures in 
the planets' atmospheres are lower than those 
in known polarized brown dwarfs and, given a 
certain temperature gradient, lower temperatures 
will yield higher P. 

A detection of an infrared polarized signal will 
immediately confirm the presence of scattering 
particles in the planetary atmosphere. In addi- 
tion, the polarization angle x will reveal the com- 
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ponents of the planet's spin axis normal to the 
line-of-sight for a zonally symmetric planet. For 
such planets, both P and F will show little vari- 
ability. On the other hand, if the polarization is 
caused by moving clouds or hot spots, both P and 
X will vary in time, and stronger than F. Hence, a 
time series of P and \ will give insight into the un- 
derlying source, and a periodic signal could even 
yield atmospheric rotation rates. However, some 
cases, like a banded and an oblate planet, may give 
very similar flux and polarization spectra and dis- 
entangling the two cases might be very difficult. 
Knowledge of the gravity on a planet will perhaps 
help in such a case, as the gravity is strongly con- 
nected to the possible oblateness. 

Together with flux measurements, polarime- 
try can also constrain particle properties, like 
their scattering albedo and size. Using polarime- 
try, one could for instance distinguish between 
absorbing iron particles and scattering silicate 
particles, which yield simila r fits to flux spectra 
([Madhusudhan et al.l l2011al) . Fig. [4] also shows 
that the polarization spectrum is much more sen- 
sitive to e.g. cloud top heights than the flux spec- 
trum. Furthermore, the broadband polarization is 
sensitive to the atmospheric temperature gradient, 
which is not easily obtained from broadband flux 
measurements. 

We thank Wiel Wauben making his scattering 
code available to us. We thank the anonymous 
referee for useful suggestions. We acknowledge fi- 
nancial support by the Netherlands Organisation 
for Scientific Research (NWO). 
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Fig. 1. — Sketch of thermal radiation with emis- 
sion angle 6, that is scattered by a particle high in 
the planetary atmosphere, with dark areas indicat- 
ing higher temperatures and thus larger thermal 
fluxes. Panel (a) illustrates a positive temperature 
gradient, and panel (b) a negative one. 
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Fig. 2. — Normalized polarized fluxes —Q/F at 
A = 1.05/xm (solid lines) and 1.11 /im (dashed 
lines) for atmospheres with (a) different tempera- 
ture gradients and a high altitude cloud with opti- 
cal thickness 0.1, and (b) a temperature gradient 
of 300 K/\np and different cloud optical thick- 
nesses. Thick lines indicate emission angles of 80° 
and thin lines those of 30°. 
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Fig. 3. — Simulated images of normalized —Q/F 
across the planetary disk for the model planets 
presented in this paper. Integrating —Q/F over 
the disk yields P — for the spherically symmetric 
planet, while it will usually leave a net P for the 
other cases. 
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Fig. 4. — Normalized F and P as a function of 
wavelength for a homogeneous, flattened planet 
with oblatcncss of 0.3. Solid lines indicate models 
with clouds at the top of the atmosphere, dotted 
lines with clouds at 0.015 bar. P is only plotted for 
F > 4 • 10~ 3 as small integration errors at these 
low flux levels give rise to very noisy P spectra 
there. 
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Fig. 5. — Normalized F, P, and \ at A 2 of a ro- 
tating planet with a 20°x20° dusty hot spot on 
its equator. The longitude denotes the distance of 
the spot with respect to the sub-observer point. 
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